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We show that the phenomenology of isospin effects on heavy ion reactions at interme- 
diate energies (few AGeV range) is extremely rich and can allow a "direct" study of the 
O ' covariant structure of the isovector interaction in the hadron medium. We work within a 



o 



relativistic transport frame, beyond a cascade picture, consistently derived from effective 



\0 • Lagrangians, where isospin effects are accounted for in the mean field and collision terms. 
^ ' Rather sensitive observables are proposed from collective flows ("differential" flows) and 
from pion/kaon production [tt^/ti^, K^/K^ yields). For the latter point relevant non- 
equilibrium effects are stressed. The possibility of the transition to a mixed hadron-quark 
phase, at high baryon and isospin density, is finally suggested. Some signatures could 
Pi . come from an expected "neutron trapping" effect. 

^ ; 1. Introduction 

\ Recently the development of new heavy ion facilities (radioactive beams) has driven 
the interest on the dynamical behaviour of asymmetric matter, see the recent reviews 
[ m 12] . Here we focus our attention on relativistic heavy ion collisions, that provide a 
unique terrestrial opportunity to probe the in-medium nuclear interaction in high density 
and high momentum regions. An effective Lagrangian approach to the hadron interacting 
system is extended to the isospin degree of freedom: within the same frame equilibrium 
properties {EoS, [ 3 ) and transport dynamics [IHE] can be consistently derived. 

Within a covariant picture of the nuclear mean field, for the description of the symmetry 
energy at saturation {a^ parameter of the Weizsaecker mass formula) (a) only the Lorentz 
vector p mesonic field, and (b) both, the vector p (repulsive) and scalar S (attractive) 
effective fields [ IHl EI can be included. In the latter case the competition between scalar 
and vector fields leads to a stiffer symmetry term at high density [IHl I2]- We present 
here observable effects, in fact enhanced, in the dynamics of heavy ion collisions. Here 
we focus our attention on collective isospin flows, in particular the elliptic ones, and on 
the isospin content of particle production, in particular kaons. We finally show that in 
the compression stage of isospin asymmetric collisions we can enter a mixed deconfined 
phase, if the EoS conditions for the existence of quark stars are met. 
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2. Relativistic Transport 

The starting point is a simple phenomenological version of the Non-Linear (with re- 
spect to the iso-scalar, Lorentz scalar a field) Walecka effective theory which corresponds 
to the Hartree or Relativistic Mean Field {RMF) approximation within the Quantum- 
Hadro-Dynamics [0]. According to this model the presence of the hadronic medium leads 
to effective masses and momenta M* = M + E^, k*'^ = — T^^, with S^, S'^ scalar 
and vector self-energies. For asymmetric matter the self-energies are different for protons 
and neutrons, depending on the isovector meson contributions. We will call the corre- 
sponding models as NLp and NLp6, respectively, and just NL the case without isovector 
interactions. For the more general NLp6 case the self-energies of protons and neutrons 
read: 

n) = -fa(T{ps) ± fspss, n) = f^f =F /pj^, {upper signs for neutrons), (1) 

where p^ = psp + psn, = jp + Ps3 = Psp - psn, js = jp - in the total and isospin 
scalar densities and currents and fa,uj,p,5 are the coupling constants of the various mesonic 
fields. (t{Ps) is the solution of the non linear equation for the a field [Elll]- 

For the description of heavy ion collisions we solve the covariant transport equation of 
the Boltzmann type [ HI E] within the Relativistic Landau Vlasov [RLV) method, using 
phase-space Gaussian test particles [^, and applying a Monte-Carlo procedure for the 
hard hadron collisions. The collision term includes elastic and inelastic processes involving 
the production/absorption of the A(1232MeV^) and A^*(1440MeK) resonances as well as 
their decays into pion channels, [p. 

It is worth to note that the nucleon mean field (Vlasov) propagation is given by the 
following equations of motion for the test particles trajectories [12]: 

A^/^ - —n^" - P^"^'^^ F^" ir■iT^'\ + d^M*ir^ (2\ 

dr""' ~ M*{x) ' dr^' ~ MtixY' l^^l^JJ + C {x) . ^Z) 

In order to have an idea of the dynamical effects of the covariant nature of the interacting 
fields, we write down, with some approximations, the "force" acting on a particle. Since 
we are interested in isospin contributions we will take into account only the isovector part 
of the interaction [ITU]: 



±/pT7I - d'^Jz Tfs^Pss ~ ±fpTh^p3Tf5^Ps3, (upper signs proton) {?,) 



dr -"^M* 



E* 



The Lorentz force (first term of Eq.®) shows a = jpr boosting of the vector coupling, 

i 

while from the second term we expect a 7-quenched 6 contribution. 

3. Collective Flows 

The fiow observables can be seen respectively as the first and second coefficients of 
a Fourier expansion of the azimuthal distribution [Ej: ^{y,Pt) ~ 1 + 2Vicos{(f)) + 

2V2Cos(20) where pt = ^pl + is the transverse momentum and y the rapidity along 
beam direction. The transverse flow can be expressed as: Viiy^pt) = (^). The sideward 
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Figure 1. Differential neutron-proton flows for the ^^'^Sn +1^4 reaction at 1.5 AGeV {b = 
6/m) from the two different models for the isovector mean fields. Left: Transverse Flows. Right: 
Elliptic Flows. Full circles and solid line: NLp5. Open circles and dashed line: NLp. 

(transverse) flow is a deflection of forwards and backwards moving particles, within the 
reaction plane. The second coefficient of the expansion defines the elliptic fiow given 

2 2 

by V2{y,pt) = ( ^^p/'^ ). It measures the competition between in-plane and out-of-plane 
emissions. The sign of V2 indicates the azimuthal anisotropy of emission: particles can 
be preferentially emitted either in the reaction plane (V2 > 0) or out-of-plane {squeeze — 
out, V2 < 0) [ im 112]. For the isospin effects the neutron-proton dif ferential flows 

^\y-,Pt) = ^12(1/ yPt) ~ ^i,2{y-iPt) have been suggested as very useful probes of the 
isovector part of the EoS since they appear rather insensitive to the isoscalar potential 
and to the in medium nuclear cross sections, [I13j. 

In heavy- ion collisions around lAGeV with radioactive beams, differential flows will 
directly exploit the Lorentz nature of a scalar and a vector fleld, see the different 7- 
boosting in the local force, Eq.®. In Fig^ transverse and elliptic differential flows are 
shown for the ^^"^Sn +^'^^ Sn reaction at 1.5 AGeV {b = Qfm), The effect of the 

different structure of the isovector channel is clear. Particularly evident is the splitting 
in the high pt region of the elliptic flow. In the (p + 6) dynamics the high-pt neutrons 
show a much larger squeeze — out. This is fully consistent with an early emission (more 
spectator shadowing) due to the larger p-fleld in the compression stage. We expect similar 
effects, even enhanced, from the measurements of differential flows for light isobars, like 
vs. ^He. 



4. Isospin effects on sub-threshold kaon production at intermediate energies 

Kaon production has been proven to be a reliable observable for the high density EoS 
in the isoscalar sector [EIEIE] Here we show that the K^'^ production (in particular 
the / yield ratio) can be also used to probe the isovector part of the EoS. 

Using our RMF transport approach we analyze pion and kaon production in central 
^^'^Au +^^^ Au collisions in the 0.8 — 1.8 AGeV beam energy range, comparing models 
giving the same "soft" EoS for symmetric matter and with different effective fleld choices 
for Esym- We will use three Lagrangians with constant nucleon-meson couplings {NL... 
type, see before) and one with density dependent couplings {DDF, see [E]), recently 
suggested for better nucleonic properties of neutron stars [Ej- In the DDF model the fp 
is exponentially decreasing with density, resulting in a rather "soft" symmetry term at high 
density. The hadron mean fleld propagation, which goes beyond the "collision cascade" 
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Figure 2. Time evolution of the A^''^'^+ resonances and pions vr^'^ (left), and kaons {K^'^ 
(right) for a central (6 = fm impact parameter) Au+Au collision at 1 AGeV incident energy. 
Transport calculation using the NL, NLp, NLp6 and DDF models for the iso-vector part of the 
nuclear EoS are shown. 

picture, is essential for particle production yields: in particular the isospin dependence of 
the self-energies directly affects the energy balance of the inelastic channels. 

Fig. 121 reports the temporal evolution of A^'^'^"*" resonances, pions (vr^'^) and kaons 
(iir+'°) for central Au+Au collisions at lAGeV. It is clear that, while the pion yield 
freezes out at times of the order of 50 fm/c, i.e. at the final stage of the reaction (and at 
low densities), kaon production occur within the very early (compression) stage, and the 
yield saturates at around 20 fm/c. From Fig. |21we see that the pion results are weakly 
dependent on the isospin part of the nuclear mean field. However, a slight increase 
(decrease) in the tt~ (tt^) multiplicity is observed when going from the NL (or DDF) to 
the NLp and then to the NLp6 model, i.e. increasing the vector contribution fp in the 
isovector channel. This trend is more pronounced for kaons, see the right panel, due to 
the high density selection of the source and the proximity to the production threshold. 

When isovector fields are included the symmetry potential energy in neutron-rich matter 
is repulsive for neutrons and attractive for protons. In a HIC this leads to a fast, pre- 
equilibrium, emission of neutrons. Such a mean field mechanism, often referred to as 
isospin fractionation [dl21, is responsible for a reduction of the neutron to proton ratio 
during the high density phase, with direct consequences on particle production in inelastic 
A^A^ collisions. 

Threshold effects represent a more subtle point. The energy conservation in a hadron 
collision in general has to be formulated in terms of the canonical momenta, i.e. for a 
reaction 1-1-2 —>■ 3-1-4 as Sin = {ki+k2)^ = (^3+^4)^ = Sout- Since hadrons are propagating 
with effective (kinetic) momenta and masses, an equivalent relation should be formulated 
starting from the effective in- medium quantities k*'^ = — H^ and m* = m + S^, where 
Ss and Tj^ are the scalar and vector self-energies, Eqs.(P3). The self-energy contributions 
will influence the particle production at the level of thresholds as well as of the phase 
space available in the final channel 

In neutron-rich colliding systems Mean field and threshold effects are acting in opposite 
directions on particle production and might compensate each other. As an example, nn 
collisions excite A~'° resonances which decay mainly to vr". In a neutron-rich matter 
the mean field effect pushes out neutrons making the matter more symmetric and thus 
decreasing the vr" yield. The threshold effect on the other hand is increasing the rate of 
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Figure 3. 238 _j_238 jj^ ^ AGeV, semicentral. Correlation between density, temperature, 
momentum thermalization inside a cubic cell, 2.5 fm wide, in the center of mass of the system. 

7r~'s due to the enhanced production of the A~ resonances: now the nn — > pA^ process 
is favored (with respect to pp nA++) since more effectively a neutron is converted into 
a proton. Such interplay between the two mechanisms cannot be fully included in a non- 
relativistic dynamics, in particular in calculations where the baryon symmetry potential 
is treated classically [UHllini- 

We have to note that in a previous study of kaon production in excited nuclear matter 
the dependence of the K^/K~^ yield ratio on the effective isovector interaction appears 
much larger (see Fig.8 of ref. [E]). The point is that in the non-equilibrium case of a 
heavy ion collision the asymmetry of the source where kaons are produced is in fact 
reduced by the n ^ p "transformation" , due to the favored nn pA" processes. This 
effect is almost absent at equilibrium due to the inverse transitions, see Fig. 3 of ref.[inj. 
Moreover in infinite nuclear matter even the fast neutron emission is not present. This 
result clearly shows that chemical equilibrium models can lead to uncorrect results when 
used for transient states of an open system. 

5. Testing Deconfinement at High Isospin Density 

The hadronic matter is expected to undergo a phase transition into a deconfined phase of 
quarks and gluons at large densities and/or high temperatures. On very general grounds, 
the transition's critical densities are expected to depend on the isospin of the system, but 
no experimental tests of this dependence have been performed so far. Moreover, up to 
now, data on the phase transition have been extracted from ultrarelativistic collisions, 
when large temperatures but low baryon densities are reached. In order to check the 
possibility of observing some precursor signals of some new physics even in collisions of 
stable nuclei at intermediate energies we have performed some event simulations for the 
collision of very heavy, neutron-rich, elements. We have chosen the reaction ^ssf/ _j_238 jj 
(average proton fraction Z/A = 0.39) at 1 AGeV and semicentral impact parameter 
b = 7 fm just to increase the neutron excess in the interacting region. 

In Fig. 01 we report the evolution of momentum distribution and baryon density in a 
space cell located in the cm. of the system. We see that after about 10 fm/c a nice 
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local equilibration is achieved. We have a unique Fermi distribution and from a simple 
fit we can evaluate the local temperature. We note that a rather exotic nuclear matter is 
formed in a transient time of the order of 10 fm/c, with baryon density around 3 — 4po, 
temperature 50 — 60 MeV, energy density 500 MeV fm~^ and proton fraction between 
0.35 and 0.40, likely inside the estimated mixed phase region, see the following.. 

Here we study the isospin dependence of the transition densities [ 120] in a systematic 
way, exploring also the possibility of forming a mixed-phase of quarks and hadrons in 
experiments at energies of the order of a few GeV per nucleon. Concerning the hadronic 
phase, we use the relativistic non-linear model of Glendenning-Moszkowski (in particular 
the "soft" GM3 choice) [ where the isovector part is treated just with p meson 
coupling, and the iso-stiffer NLpS interaction For the quark phase we consider the 

MIT bag model [l23j with various bag pressure constants. In particular we are interested 
in those parameter sets which would allow the existence of quark stars [ |^ 123; i-e. 
parameters sets for which the so-called Witten-Bodmer hypothesis is satisfied [ EBl I27j. 
One of the aim of our work it to show that if quark stars are indeed possible, it is then 
very likely to find signals of the formation of a mixed quark-hadron phase in intermediate- 
energy heavy- ion experiments [I22j. 

The structure of the mixed phase is obtained by imposing the Gibbs conditions [EDI EI] 
for chemical potentials and pressure and by requiring the conservation of the total baryon 
and isospin densities 

PB = {1- x)Pb + XpI , P3 = (1 - X)P3 + XP? , (4) 

where x is the fraction of quark matter in the mixed phase. In this way we get the hinodal 
surface which gives the phase coexistence region in the (T,pB,P3) space [EH 120]- For a 
fixed value of the conserved charge ps we will study the boundaries of the mixed phase 
region in the (T, p^) plane. In the hadronic phase the charge chemical potential is given 
by = ^Esyjnlps)— ■ Thus, we expect critical densities rather sensitive to the isovector 
channel in the hadronic EoS. 

In Fig. m we show the crossing density pcr separating nuclear matter from the quark- 
nucleon mixed phase, as a function of the proton fraction Z/A. We can see the effect of 
the (5-coupling towards an earlier crossing due to the larger symmetry repulsion at high 
baryon densities. In the same figure we report the paths in the (p, Z/A) plane followed 
in the cm. region during the collision of the n-rich ^^^Sn-|-^'^^Sn system, at different 
energies. At 300 AMeV we are just reaching the border of the mixed phase, and we are 
well inside it at 1 AGeV. Statistical fluctuations could help in reducing the density at 
which drops of quark matter form. The reason is that a small bubble can be energetically 
favored if it contains quarks whose Z/A ratio is smaller than the average value of the 
surrounding region. This is again due to the strong Z/A dependence of the free energy, 
which favors clusters having a small electric charge. Moreover, since fluctuations favor the 
formation of bubbles having a smaller Z/A, neutron emission from the central collision 
area should be suppressed, which could give origin to specific signatures of the mechanism 
described in this paper. This corresponds to a neutron trapping effect, supported also by a 
symmetry energy difference in the two phases. In fact while in the hadron phase we have 
a large neutron potential repulsion (in particular in the NLp6 case), in the quark phase 
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Figure 4. Variation of the transition density with proton fraction for various hadronic EoS 
parameterizations. Dotted hne: GM3 parametrization; dashed Une: NLp parametrization; soUd 
Une: NLp6 parametrization. For the quark EoS, the MIT bag model with 5^/^=150 MeV. 
The points represent the path fohowed in the interaction zone during a semi-central ^^^Sn+^^^Sn 
collision at 1 AGeV (circles) and at 300 AMeV (crosses). 

we only have the much smaller kinetic contribution. If in a pure hadronic phase neutrons 
are quickly emitted or "transformed" in protons by inelastic collisions, when the mixed 
phase starts forming, neutrons are kept in the interacting system up to the subsequent 
hadronization in the expansion stage • Observables related to such neutron "trapping" 
could be an inversion in the trend of the formation of neutron rich fragments and/or of 
the 7r~/7r"'", K^/K~^ yield ratios for reaction products coming from high density regions, 
i.e. with large transverse momenta. In general we would expect a modification of the 
rapidity distribution of the emitted "isospin", with an enhancement at mid-rapidity joint 
to large event by event fluctuations.. 

6. Perspectives 

We have shown that collisions of n-rich heavy ions at intermediate energies can bring 
new information on the isovector part of the in-medium interaction at high baryon den- 
sities, qualitatively different from equilibrium EoS properties. We have presented quan- 
titative results for differential collective flows and yields of charged pion and kaon ratios. 
Important non-equilibrium effects for particle production are stressed. Finally our study 
supports the possibility of observing precursor signals of the phase transition to a mixed 
hadron-quark matter at high baryon density in the collision, central or semi-central, of 
neutron-rich heavy ions in the energy range of a few GeV per nucleon. As signatures 
we suggest to look at observables particularly sensitive to the expected different isospin 
content of the two phases, which leads to a neutron trapping in the quark clusters. The 
isospin structure of hadrons produced at high transverse momentum should be a good 
indicator of the effect. 

In conclusion the results presented here appear very promising for the possibility of ex- 
tracting information from terrestrial laboratories on the Lorentz structure of the isovector 
nuclear interaction in a medium at densities of astrophysical interest. The use of radioac- 
tive beams at relativistic energies would be extremely important. 
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hadron-quark phase transition at high baryon and isospin density. 
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